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Abstract. This paper presents a detailed study of the chemistry in the outflow associated with the low-mass 
protostar NGC 1333-IRAS2A down to 3" (650 AU) scales. Millimeter-wavelength aperture-synthesis observations 
from the Owens Valley and Berkeley-IUinois-Maryland- Association interferometers and (sub)millimeter single- 
dish observations from the Onsala Space Observatory 20 m telescope and Caltech Submillimeter Observatory are 
presented. The interaction of the highly collimated protostellar outflow with a molecular condensation ~15000 AU 
from the central protostar is clearly traced by molecular species such as HCN, SiO, SO, CS, and CH3OH. Especially 
SiO traces a narrow high velocity component at the interface between the outflow and the molecular condensation. 
Multi-transition single-dish observations are used to distinguish the chemistry of the shock from that of the 
molecular condensation and to address the physical conditions therein. Statistical equilibrium calculations reveal 
temperatures of 20 and 70 K for the quiescent and shocked components, respectively, and densities near lO*" cm~^. 
The line-profiles of low- and high-excitation lines are remarkably similar, indicating that the physical properties 
are quite homogeneous within each component. Significant abundance enhancements of two to four orders of 
magnitude are found in the shocked region for molecules such as CH3OH, SiO and the sulfur-bearing molecules. 
HCO"^ is seen only in the aftermath of the shock consistent with models where it is destroyed through release 
of H2O from grain mantles in the shock. N2H''' shows narrow lines, not affected by the outflow but rather 
probing the ambient cloud. The overall molecular inventory is compared to other outflow regions and protostellar 
environments. Differences in abundances of HCN, H2CO and CS are seen between different outflow regions and 
are suggested to be related to differences in the atomic carbon abundance. Compared to the warm inner parts of 
protostellar envelopes, higher abundances of in particular CHaOIf and SiO are found in the outflows, which may 
be related to density differences between the regions. 
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1. Introduction 

One of the manifestations of a newly formed low-mass 
protostar is the presence of a highly collimated and ener- 
getic outflow or jet. A natural consequence of the propa- 
gation of such high velocity outflows through the proto- 
stellar envel ope and the ambient mo lecular medium are 
shockfronts ijReipurth fc Raeal Il999|) . Shocks both heat 
and compress the gas and also trigger chemical reactions 
in the gas-phase, leading to a different chemistry than 
observed otherwise. Shock processing of dust grains may 
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lead to the injection of atoms and molecules back into 
the gas, which further distinguishes the chemistry in the 
shocked region from that of a quiescent protostellar envi- 
ronment fe.g. . iTielenSi.1999) . This paper presents a study 
of the physics and chemistry of the outflow associated 
with a well-known class young stellar object NGC 1333- 
IRAS2 using high-resolution millimeter wavelength inter- 
ferometer and multi-transition single-dish observations. 
The NGC 1333-IRAS2 outflow provides a unique oppor- 
tunity to study the effects of outflows on ambient molec- 
ular clouds, as the main shock is well separated from the 
central protostar and shows a relatively simple morphol- 
ogy. The combination of single-dish and interferometry 
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observations makes it possible to discuss the physical and 
chemical properties of the outflowing gas and to address 
the spatial differentiation of the chemistry in the outflow 
region resolved by the interferometer observations. 

Studies of molecular abundances in regions of high 
outflow activity provide insight into the dependence of 
the chemical reaction networks on temperature and den- 
sity. Furthermore it is important to recognize the effect of 
outflow-triggered chemistry in the inner protostellar enve- 
lope, to disentangle it from emission from a circumstellar 
disk or to address the effect of passive heating by the cen- 
tral protostar. In the central part of the protostellar enve- 
lope, thermal evaporation of dust grain mantles can lead 



protostars (e.g IRAS 16293-2422, 


Ceccarelh et al.l2000albl 


ISchoier et alJl2002l ICazaux et al. 


l2nn,'*. 



NGC 1333-IRAS2 (also known as IRAS 03258-f3104; 
hereafter simply IRAS2) is located in the NGC 1333 
molecular cloud, harboring several class and I objects, 
first i dentified through IRAS maps by l.lennings et alJ 
(^13) ■ Continuum observations reveal that IRAS2 is a 
binary source with two com ponents, IRAS2A and 2B, 
separ a ted bv 6500 AU (3 0'^) llSandell et al.lll994t lBlak3 
Il996l iLoonev et all 1200(1) . IRAS2A is responsible for 
a highly coUimated east-west outflow giving rise to a 
strong shock ~ 15000 AU from the central continuum 
source (Fig. A strong CO outfl ow in the north- 
south direction has also been o bserved ijLiseau et alll988t 
lEneargiola &: Plambecklll999(), which originates with in a 
few arcseconds from IRAS2A l|,Te)rgensen et al.ll2n03aj) . 

iLanger et alJ l)l996|) mapped the entire NGC 1333 re- 
gion in CS and identifled two peaks in CS emission toward 
the IRAS2 outflow. They suggested that these are asso- 
ciated with red-shifted (eastern) and blue- shifted (west- 
ern) b ow shock components of the outflow. ISandell et al.l 
lll994^ reported bright CII3OH emission toward the red- 
shifted (eastern) outflow lobe and correspondingly a high 
CH3OII abundance, possibly enhanced by the shock. 
The more detailed structure of the CH3OH emission 
from the IRAS2A e ast-west outflow was discussed by 
iBachiller et al.l l)l998j) . who mapped the outflow positions 
at 3" using the IRAM inter ferometer and 30 m single- 
dish telescope. IBachiller et al.l derived the physical condi- 
tions in the shock interaction zone from LVG calculations 
and obtained a density of ^ 10^ cm^^ and temperature of 
~ 100 K. 

IBachiller et alJ also found that the observed methanol 
emission translates to a large enhancement of CII3OII by 
a factor 300 in the IRAS2A outflow. CII3OII is thought 
to be released directly from the dust grain mantles and 
is oft en seen to be associa ted with protostellar outflows 
fe.g.. ISachiller et al.lll995l) . Other often- used tracers of 
shocks associated with protostellar outflows are Si-bearing: 
species, in particular S i O (e.g., iMartin-Pintado et alJ 



around low-mass protostars llBac hiller &: Perez Gutierred 
ll997HGarav et alJll998t fBachiller et alJl200lh . 

In this paper we present a study of the detailed chem- 
istry of the shock associated with the IRAS2A outflow 
based on observations of a wide range of molecular lines 
at ~ 3 — 6" resolution from the Owens Valley Radio 
Observatory (OVRO) and Berkeley Maryland Illinois 
Association (BIMA) millimeter interferometers, together 
with millimeter and submillimeter single dish observations 
from the Onsala 20 m telescope (OSO) and the Caltech 
Submillimeter Observatory 10.4 m telescope (CSO). Parts 
of the OVRO o bservations have previously been presented 
bv lBlakel l)l996|) . Sect. |2| describes the observations and re- 
ductions. The maps from the interferometry observations 
are presented and discussed in Sect. 13.11 while the single 
dish observations are treated in Sect. 13.21 The physical 
and chemical properties of the shock region are analyzed 
using statistical equilibrium calculations as described in 
Sect.0]and molecular abundances are derived. Sect.jSjdis- 
cusses the inferred chemistry and compares it to other 
well-studied outflow regions, to other types of star- forming 
environments and to available models for the chemistry 
in outflow regions. The main fl ndings are summarized 
in Sect. |B| A companion paper ijjorgensen et al.ll2003a(l 
presents details of a millimeter-wavelength interferometer 
study of the environment surrounding the central proto- 
stellar system. 



2. Overview of observations 

The position of the shock in the eastern lobe of the out- 
flow associated with IRAS2A (a(2000) = 03^29™0ffi0; 
(5(2000) = 31°14'19'.'0) was observed with the Milhmeter 
Array of the Owens Valley Radio Observatory (OVRO)"'^ 
between October 5, 1994 and January 1, 1995 in the six- 
antenna L- and H-configurations. Tracks were obtained in 
two frequency settings at 86 and 97 GHz, and each track 
observed in alternatingly two fields: the bow shock at the 
end of the eastern outflow discussed in th is paper and the 
positio n of the central protostellar source llorgen sen et alJ 
l2003al) . The observed tracks cover projected baselines of 
3.1-70 kA at 86 GHz and the lines observed are listed in 
Table ^ The lines were recorded in spectral bands with 
widths of 32 MHz (- 100 km s^^). }1^^G0+ 1 - and CS 
2 — 1 were observed in 128 spectral channels, the remain- 
ing line setups included 64 spectral channels. The complex 
gain variations were calibrated by observing the nearby 
quasars 0234-f 285 and 3C84 approximately every 20 min- 
utes. Fluxes were calibrated by observations of Uranus 
and Neptune. The rms noise levels are 0.05 Jy beam"-'^ 
in the 250 kHz channels with a synthesized beam size 
of 3.2 X 2.8". Calibration and flagging of visibilities with 
clearly deviating amplitudes and/o r phases was perform ed 
with the MMA reduction package ijScoville et al.lll99^ . 



ll992HCodella et al.lll999HGarav et al.ll2002^ . High abun- 
dances of these species may mark a clear distinction 
of shocked gas from unprocessed gas in the envelopes 



^ The Owens Valley Millimeter Array is operated by the 
California Institute of Technology under funding from the US 
National Science Foundation (grant no. AST-9981546). 




Fig. 1. Overview of the IRAS2A outflow region. The colored image shows the SCUBA 850 /im emission tracing the 
cold dust. N2H+ BIMA line observations and SiO O VRO line observations (this pa per) are indicated by blue and 
green contours, respectively. CO 2-1 emission from the'Enga rgiola fc: Plambeca l)l99 9^ is indicated by the black/white 
dashed contours. The red stars indicate the position of IRAS2A and IRAS2B from iJorgensen ct al, (,2003a) . 



The millimeter interferometer of the Berkeley-Illinois- 
Maryland Association (BIMA)^ observed the IRAS2A 
outflow position between March 4 and April 15, 2003. 
The array B- and C-configurations provided projected 
baselines of 2.7-71 kA. The hues of HCO+ 1-0, HCN 
1-0, N2H+ 1-0, and C^^S 2-1 were recorded in 256- 
channel spectral bands with a total width of 6.25 MHz 
(r^ 20 km s~^). The complex gain of the interferometer was 
calibrated by observing the bright quasars 3C84 (4.2 Jy) 
and 0237-1-288 (2.3 Jy) approximately every 20 minutes. 
The absolute flux scale was bootstrapped from observa- 
tions of Uranus. The rms noise levels are 0.2 Jy beani^^ 
in the 24 kHz channels, with a synthesized beam size of 



^ The BIMA array is operated by the Universities of 
California (Berkeley), Illinois, and Maryland, with support 
from the National Science Foundation (grants AST-9981308, 
AST-9981363 and AST-9981289). 



6.1" X 5.0" FWHM (7.6" x 6.8" for the C^^S and N2H+ ob- 
servations). The data were calibrated with routines from 
the MIRIAD software package (Sault et al. 1995). 

In addition to the interferometry data, a num- 
ber of molecular lines were observed toward the po- 
sition of the red-shifted shock using the Caltech 
Submillimeter Observatory 10.4 m (CSO)^ and Onsala 
Space Observatory 20 m (OSO)** telescopes. The point- 
ing was checked regularly and found to be accurate to a 
few arcseconds. The typical beam sizes are 45"-33" for the 
OSO 20 m (86-115 GHz) and 26"-20" for the CSO (217- 



^ The Caltech Submillimeter Observatory 10.4 m is oper- 
ated by Caltech under a contract from the National Science 
Foundation (grant no. AST-9980846). 

* The Onsala 20 m telescope is operated by the Swedish 
National Facility for Radio Astronomy, Onsala Space 
Observatory at Chalmers University of Technology. 
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Table 1. Overview of the observations of the IRAS2A 
outflow position treated in this paper. 



Line Rest freq. Observed with 



CO 


2 - 


1 


230.5380 


cso 




3- 


2 


345.7960 


cso 


c^^o 


3- 


2 


329.3305 


cso 


CH30H 


2i - 


- li 


97.5828 


OVRO 




72 - 


-61 


338.7222 


CSO 


cs 


2 - 


1 


97.9810 


OSO, OVRO 




5- 


4 


244.9356 


CSO 




7 - 


6 


342.8830 


cso 




2 - 


1 


96.4129 


BIMA 


HON 


1 - 


0" 


88.6318 


OSO, BIMA 




4 - 


3 


354.5055 


CSO 


HCO+ 


1 - 





89.1885 


OSO, BIMA 




4~ 


3 


356.7343 


CSO 


H"CO+ 


1 - 





86.7543 


OSO, OVRO 


H2CO 


5l,5 


- 4i,4 


351.7686 


CSO 




5o,5 


- 4o,4 


362.7359 


CSO 


N2H+ 


1 - 


0" 


93.1737 


OSO, BIMA 


SiO 


2 - 


1 


86.8470 


OSO, OVRO 




5- 


4 


217.1049 


CSO 




8- 


7 


347.3306 


CSO 


so 


22 ' 


- li 


86.0940 


OVRO 




23 - 


- I2 


99.2999 


OSO 




89 ' 


-78 


346.5285 


CSO 


S02 


73,5 


— 82,6 


97.7024 


OVRO 



"Hyperfine splitting - multiple lines observed in one setting. 
The coordinates for the single-dish observations are Qf(2000) = 
03''29™01':0, 5(2000) = 31°14'20", i.e., corresponding to offset 
of (79",-17") in the maps presented in this paper. 

356 GHz) observations. The data were calibrated using 
the standard chopper wheel method. The spectra were re- 
duced in a standard way by subtracting baselines and by 
dividing by the main-beam efficiencies Tymb as given on the 
web pages for the two telescopes, r^mb ranges from 0.6 to 
0.43 for frequencies of 86 to 115 GHz for the OSO 20 m 
and 0.67 to 0.62 for frequencies of 217 to 356 GHz for the 
CSO. An overview of all the observed lines (single-dish 
and interferometer) is given in Tabled 

3. Data 

3.1. Interferometry 

Fig. I2I3I show moment maps for the lines observed at 
OVRO (CS, SO, SiO and CH3OH) and BIMA (i.e., 
HCO+, HCN, N2H+ and C^^S). In all maps the coor- 
dinates are given as offsets relative to the position of 
the central protostar, IRAS2A: a(2000) = 03'^28'"55!7, 
(5(2000) = 31°14'37". Emission of SO2 and Hi3C0+ was 
not detected in the interferometer maps toward the out- 
flow position. 

Most of the observed lines show indications of material 
affected by the outflow with clear line wings spreading out 
to 10-15 km s^^ from the systemic velocity of 7 km s~^. 



One exception is N2H'^ which shows narrow hyperfine 
components of approximately 1 km s^^ width (FWHM). 
Of the observed lines, SiO, CS and HCN show emis- 
sion stretching furthest from the systemic velocity (out to 
« 20 km s^i) while the remaining species show somewhat 
narrower profiles (wings stretching out to « 10 km s~^ rel- 
ative to the systemic velocity), with HCO"*" showing most 
material closest to the cloud systemic velocity. SiO is not 
seen at low velocities as is also the case for CH3OH and 
SO. 

In general the molecules observed at BIMA most 
clearly trace the extended low velocity material. This may 
in part be due to the different (u, v) coverage of the two 
arrays, with the BIMA observations recovering more of 
the extended emission, but the species observed at OVRO 
may also be those that are predominantly enhanced by 
the outflow generated shock. 

The N2H+ emission traces a ridge of material with a 
number of "cores" stretching from the north-east of the 
map toward the center and back again to the north-west. 
A dominant core is seen close to the center (offsets of 
(67", -3")) which is also picked up by the HCO"'' maps. 
South of this core a "<" -shaped extension is seen in both 
HCO"'" and N2H+. The low velocity CS emission picks up 
only this feature. A similar component was also s een in 
the CH3OH emission mapped by Bachiller et alJ l(l998^ 
but is not evident in the CH3OH observations presented 
here, possibly due to lower sensitivity. 

The high-velocity material is generally much less ex- 
tended than the low- velocity material. The HCN, SiO and 
CS trace a narrow component stretching 30-40" along the 
outflow propagation and 5-10" in the direction perpendic- 
ular to this. The narrow component points directly to the 
"<" -shaped feature in the low velocity emission material. 
The HCN emission is slightly more extended than that 
of the two other species, again likely due to the different 
(it, v) coverage of the observations from the two arrays. 
SO and CH3OH show a slightly weaker structure along 
the same narrow component. 

The emission of CH3OH and SO is located downstream 
(west) of the outflow propagation direction compared to, 
e.g., the peak of SiO. Even farther downstream around 
offset (59",-13"), HCN, CS and SiO show another strong 
feature where the N2H"'" emission "pinches" the outflow. 
The HCO"'" wing emission is seen only at this position and 
is found to be more extended, filling out the region void 
of N2H"'". In fact the HCO"'" emission can be traced all the 
way back to the central protostar as is also the case for 
CO (Fig.lIJ. It is striking how the HCO"^ and, e.g., HCN 
wing emission trace significantly different components, im- 
plying a clear chemical differentiation. 

Position-velocity diagrams for CS and SiO are pre- 
sented in Fig. |S1 Note the symmetry around the X-axis 
in these diagrams with low- velocity emission constituting 
a broad component of weak emission. For both species, 
the high velocity component is more pronounced toward 
the working surface of the outflow. 
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Fig. 2. Moment maps for the lines observed at OVRO: a) CS, b) SO, c) SiO and d) CH3OH. The black lines indicate 
low velocity emission integrated over velocities of 5-9 km s^^ whereas the red lines indicate high velocity emission 
integrated over 9-16 km s~^ for CH3OH and SO and 9-25 km s~^ for CS and SiO. The contours are given in steps 
of 3(7 and are overlaid on grey-scale images of the low velocity N2H"*' emission. The x- and y-axis offsets are given 
relative to the position of the IRAS2A source. In the upper left panel the arrow indicates the direction back to the 
central protostar. 



3.2. Single-dish 

The observed single-dish spectra are presented in Fig. 
In agreement with the interferometry maps, the N2H"'" hy- 
perfine lines show Gaussian profiles with no sign of outflow 
wings. The same applies for the C^'^O 3-2 observations. 
The 3 mm lines of HCO+, HON, CS and SO show a "two 
component" line profile with a narrow peaked profile close 
to the systemic velocity of the cloud along with a clear 



red wing extension. SiO on the other hand does not show 
the narrow component but has a more or less abrupt in- 
crease slightly above the cloud rest velocity with a close to 
linear decline in strength toward the higher (red-shifted) 
velocities. For the CSO (0.8-1.4 mm) lines a similar trend 
is seen (except for C^^O 3-2). For the CO 2-1 and 3-2 
transitions two strong dips are seen at the cloud rest ve- 
locity and in the blue part of the spectrum likely due to 
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Fig. 3. Moment maps for the lines observed at BIMA: a) HCN, b) HCO+, c) N2H+ and d) C^^S. As in Fig. the 
black and red lines indicate low and high velocity material, respectively (the low velocity emission is integrated over 
5-9 km and the high velocity emission over 9-16 km s~^). The contours ascend in steps of 3a and are overlaid on 
grey-scale images of the low velocity N2H+ emission. 



self-absorption and possible switching onto non-source CO 
emission. This, however, does not affect the wing emission 
in the spectrum, which stretches out to ~ 30 km s^^ (20- 
25 km s^^ relative to the cloud rest velocity). 

In Fig. 13 spectra for individual molecules are com- 
pared. The two component separation into a core and wing 
profile seems to be unique for each molecule: the widths of 
the core part of the lines (where observable in more than 
one transition) and the dependence of line strength with 
velocity in the wings, in particular the terminal velocities 



inferred for each molecule, are both independent of the 
observed transition. This gives a clear indication that two 
distinct components with different excitation conditions 
and chemical properties are observed and that the differ- 
ent molecules probe distinct parts of each of these compo- 
nents, within which the excitation properties do not vary 
significantly. 

As for the interferometry data, the terminal veloci- 
ties seem to indicate different dynamical components of 
the outflow region. SiO and CO probe material at the 




Fig. 4. Channel maps for a) HCN, b) HCO+, c) CS and d) SiO. Contours are given in steps of 3a. The synthesized 
beam is indicated in the bottom left panel of each figure. 



highest velocities relative to the cloud rest velocity of 
20-25 km s~^, whereas HCO'*' has the lowest degree of 
wing emission extending < 10 km s^^ from the rest ve- 
locity. The remaining molecules fall somewhere in be- 
tw een. A similar trend was seen in the L1157 outflow 
bv lBachiller et alJ l)200l|) . who also found CO and SiO to 
have significantly higher terminal velocities than HCO"*", 
H2CO, SO and CS. They suggested that this could be 
related to differenc es in the forma tion mechanisms for 
the various species. iBachiller et al .^ also argued in favor 
of rather homogeneous excitation conditions in the L1157 
outflow, since maps of various transitions for specific 



molecules were found to be very similar. In this case the 
main differences among different molecular species would 
be the result of the chemistry in the region. 

3.3. Qualitative scenario 

The morphology of the interferometry maps and the line- 
profiles of the single-dish observations can be explained in 
a simple picture as signatures of a highly coUimated out- 
flow ramming into a quiescent core or static cloud traced 
by the clumps of C^^O and N2H+ emission. The action 
of the outflow leads to sputtering of silicon off dust grains 
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Fig. 5. CS (upper) and SiO (lower) position-velocity maps for the IRAS2A shock. The coordinate frame has been 
rotated and translated to an XY-coordinate system with the X-axis along the propagation direction of the outflow 
(~ 19° with the RA axis in Fig. and the Y-axis in the perpendicular/north direction. The (0,0) point for the XY 
coordinate system has been chosen to be at (87",-23") i.e., at the working surface, or head, of the outflow as judged 
from the morphology of the high velocity emission. The contours are given at 2cr, 4(t, 8a, . . . and upwards in steps of 

4f7. 



to form SiO, which is totally absent in the quiescent core. 
At the same time the abundances of CH3OH, CS and SO 
increase, more likely due to evaporation of ice mantle ma- 
terial. A sequence in chemistry can be seen with SiO prob- 
ing material at the highest velocities followed by CS and 
SO, thought to be a result of enhanced sulfur gas phase 
chemistry, and CH3OH, H2CO and HCN, which are likely 
to be direct results of grain mantle release. The differ- 
ence in extent of the line wings indicates either that the 
molecule formation time scales are varying, with SiO be- 
ing produced most rapidly in the shock, or that the more 
volatile species such as CH3OH do not survive at the high- 
est velocities in the flow. 

That N2H+ is observed only in the quiescent cloud 
material is explained if the temperature in the material 
affected by the outflow increases to > 20 K. At this 
temperature CO is released from grain mantles and be- 
comes the dominant destruction channel of N2H"*', lower- 
ing the abundance of this molecule. Comparison between 



the morphology of the N2H+ emission and that of other 
species indicates a clear interaction between the outflow- 
ing material and the ambient cloud. Fig. I2I3I show a cav- 
ity of N2H+ emission where the shocked gas (e.g., SiO) 
appears. The two peaks seen in SiO, HCN and CS also 
seem to be related to an increase in N2H+ emission. A 
natural question is whether the outflow shapes or is be- 
ing shaped by the ambient material. Judging from the 
morphology of the larger scale emission in Fig. ^ it is 
striking to note the presence of large amounts of dust 
northeast of the central IRAS2A protostar. The two CO 
outflows identified by different authors (e.g..lLiseau et alJ 
llQSStlEngargiola fc Plambeckll99<AlKnee & Sandell200fl) 
trace the edges of this dust condensation. The two per- 
pendicular outflows seen toward IRAS2 may therefore re- 
flect the conditions in the ambient cloud material rather 
than the intrinsic properties of the central protostellar sys- 
tem: the CO outflow could simply be deflected around 
the dense material traced by the N2H+ and continuum 
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Fig. 6. Observed single dish spectra toward the IRAS2A shock position. In the upper two rows the 3 mm observations 
from the Onsala 20 m are presented, whereas the 0.8-1.4 mm observations from the CSO are shown in the lower 3 
rows. All spectra are on the Tmb scale - note that some spectra have been scaled as indicated to fit on the composite 
plot. 



emission leading to the quadrupolar morphology. In ei- 
ther case, however, this does not change the interpreta- 
tion of the shocked material in this paper. High veloc- 
ity gas (> 20 km s~^ relative to the systemic velocity) is 
present toward the eastern lobe as it is seen from most of 
the species observed in this paper and this indicates the 
presence of the shock. 

4. Analysis 

4.1. Line intensities 

Despite the clear separation between the "quiescent" and 
"shocked" parts of the line profiles for most molecules, 
the disentanglement of the emission into a core and wing 
component is not unique. Each line was decomposed in two 
parts, with the emission integrated over velocities higher 
and lower (wing and core components, respectively) than 
where the profiles in Fig.[7|separate. For C^^O and N2II+ 
Gaussians were fitted to each line giving widths of 1.5- 
2 km s~^ (FWHM), which compares well to the widths 



of the core parts of the remaining lines. Table |21 lists the 
resulting line intensities. 

The hyperfine splitting of the HCN 1-0 line gives rise 
to three components within the same setting. The two 
weaker transitions are offset —5 and 7 km s~^ relative to 
the main hyperfine line. Overlap in the line wings there- 
fore makes the interpretation of this line difficult. If the 
emitting material is in local thermodynamical equilibrium 
(LTE) and the emission is optically thin, one should ex- 
pect the hyperfine components to be in a ratio of 1:3:5. In 
Fig.|Sla comparison between the HCO+ and HCN spectra 
toward the shock position is shown. The HCN spectrum 
has been overplotted with a composite of three versions 
of the HCO+ spectrum shifted according to the frequency 
shifts of the hyperfine lines and scaled in the relative 1:3:5 
proportions - leaving an overall "normalization factor" be- 
tween the HCN and HCO+ line intensities as the only free 
parameter. The good agreement is remarkable and indi- 
cates that HCO+ and HCN trace the same material, espe- 
cially at the lower velocities. The intensity of the main hy- 
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Fig. 7. Comparison between profiles for different transitions of specific molecules. Note that the spectra have been 
smoothed by up to a factor of 10 to bring out the agreement between the lines, except for the scale factor given in the 
upper right corner of each plot. 



perfine component of the HCN 1-0 line is thereby found to 4.2. Tying interferometry and single-dish observations 
be 0.5 times the intensity of the HCO+ 1-0 line. It should together 
still be re-emphasized, that the interferometer maps show 
clear differences for the less extended HCN and HCO+ 

One question to address is whether the single-dish and 
interferometry observations trace the same material. The 



emission at higher velocities. 
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Table 2. Line parameters from single dish observations 



Molecule 


Transition 


V'(core)'' 
km s~^ 


/(core)" 
K km s-i 


y(wing)'' 
km s^^ 


/(wing)" 
K km s-i 




2-1 


[4,111 


45 


[11,30] 


13 




3-2 


[4,111 


50 


11,30 


24 




3-2 


Gaussian 


1.7 






CH3OH 


72—62 


[6,91 


2.0 


[9,211 

L ' J 


6.4 


cs 


2-1 


[5.5,9.51 


12 


[9.5,16] 


3.9 




7-6 


[5.5,9.5] 


1.6 


[9.5,16] 


1.6 


HON 


l-O'' 


[5.5,7.5] 


4.2 


[7.5,15] 


3.3 




4-3 


[5.5,7.5] 


0.72 


[7.5,15] 


0.78 


HCO+ 


1-0 


[5.5,7.5] 


8.4 


[7.5,15] 


6.5 




4-3 


[5.5,7.5] 


2.7 


[7.5,15] 


2.0 


H2CO 


5l,5 - 4l,4 


[5.5,7.5] 


3.1 


[7.5,15] 


1.5 


SiO 


2-1 






[7,29] 


10.4 




5-4 






[7.5,20] 


6.5 




8-7 






[7.5,28] 


4.4 


so 


23 - I2 


[5,9] 


4.3 


[9,17] 


4.2 




8g — 78 


[5,9] 


0.61 


[9,17] 


2.0 



Notes: "Velocity interval over which the emission is integrated, 'integrated line intensity {I = J Tmbdv). '^Main hyperfine line; 
line intensities derived through decomposition with HCO^ line. See description in text. 
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Fig. 8. The HCN 1 — spectrum toward the shock po- 
sition (black) overplotted with a composite of three ver- 
sions of the HCO"'" spectrum toward the same position 
(red) - each shifted with the measured shifts between the 
HCN hyperfine lines and scaled according (relatively) to 
the 1:3:5 line-ratio expected for optically thin emission 
from material in LTE. 

similar trends seen in the data from the two types of obser- 
vations seem to support this, but it is also known that the 
interferometry observations lack sensitivity to extended 
structures on larger scales. 

Fig. El presents comparisons between the single-dish 
observations and spectra taken from interferometry data 
sets, restored with beam sizes appropriate for the Onsala 
20 m telescope. The interferometry spectra have been 
scaled to match the wings of the single-dish spectra. Only 
a small correction of the order of 30% needs to be ap- 
plied to match the CS and SiO spectra, which is close to 
the calibration uncertainty. This confirms that the higher 



velocity emission is relatively compact (as the maps also 
suggest for the emission in the transverse direction of the 
outflow), making it less subject to the incomplete {u,v) 
sampling of the interferometer observations. 

Closer to the systemic velocity of the cloud (w 
7 km s~^) the discrepancy between the single-dish and 
interferometry spectra increases. The CS interferometry 
observations pick up only a small fraction of the emission 
in the "core" part of the single-dish line. The dip seen in 
the single-dish CS spectra at the rest velocity of the cloud 
is a result of self-absorption, while for the interferometry 
observations it is caused by the interferometer resolving 
out extended emission close to the cloud systemic veloc- 
ity. The (m, v) sampling is also responsible for the lack of 
emission in the SiO interferometry spectra at velocities 
close to the rest velocity, although it is less significant for 
this molecule. 

For the BIMA observations of HCN and HCO+, emis- 
sion close to the systemic velocity is still mostly resolved 
out as indicated by the dips in the interferometer HCN and 
HCO"*" spectra and as seen in the channel maps in Fig. 01 
The slightly better (u, v) coverage from BIMA, however, 
makes these lines less subject to resolving out at velocities 
different from the systemic velocity. 

Nyquist sampled single-dish maps of the different 
molecular species would make it possible to combine the 
interferometer and single-dish data to create maps includ- 
ing the short-spacings. This would settle the issue of the 
differing (it, v) coverage of the two arrays. In this paper 
we only have single pointing observations. The agreement 
between the single-dish and interferometry spectra in the 
line wings, however, justifies the discussion of the outflow 
component based on the morphology in the interferometry 
maps presented in Sect. 13.11 The agreement also makes 
it possible to use the interferometer maps to determine 
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Fig. 9. Comparison between the single-dish observations 
(black) and corresponding spectra from the interferome- 
ter observations restored with the single-dish beam (red) 
at the shock position. The spectra from the interferome- 
try observations have been scaled to resemble the wing of 
the single-dish spectra with the factors indicated in the 
upper right corner. For C'^^S the single-dish CS has been 
downscaled by a factor 22. 



the spatial extent of the wing component, and thereby to 
estimate the beam filling factor for the single-dish obser- 
vations. The wing part of the SiO and CS interferometry 
maps give a rough estimate of the extent of the outflow 
emission in the transverse direction of 5-10", leading to 
filling factors ranging from 0.07 to 0.32 for the single-dish 
beam sizes. For the core component of the lines on the 
other hand, the interferometer observations are less use- 
ful because of the significant fraction of the low-velocity, 
extended emission that is resolved out. For the following 
discussion, a filling factor of unity is therefore assumed for 
the core component of the single-dish data, which seems 
realistic as the interferometry maps do reveal emission ex- 
tended over scales larger than 20-30". 

4.3. Statistical equilibrium calculations 

In order to derive the physical properties of the emitting 
gas and the column densities of the various molecules, 
one-dimensional statistical equilibrium calculations were 
performed using the ID radiative transfer code called 
Radex. Radex uses the escape probability formalism to 



solve the statistical equilibrium equations f or a medium 
with constant density ar id temperature (see Ijansen et alJ 
h994ISrhmer et a,lJl2003l^ . For subthermally excited emis- 
sion this approach is an improvement compared with the 
rotation diagram method, as demonstrated for interpre- 
tations of_tlie_CH^OH emi^ toward the IRAS2A out- 
flow bv lBachiller et af] l)l998j) . The statistical equilibrium 
calculations also treats opacity effects in the correct way, 
again contrasting the rotation diagram analysis which re- 
lies on optically thin emission. 

For each molecule, line intensities were calculated for 
varying column density, density of the main collision part- 
ner (H2) and kinetic temperature. The filling factors esti- 
mated on the basis of the interferometry observations as 
discussed above were adopted and the calculated line in- 
tensities were compared to the observed ones. Since our 
main interest is in the relative behavior of the lines some 
of the uncertainties in the assumptions, e.g., the filling fac- 
tor of unity for the core component, will cancel out, if the 
differences between the observed lines are not too large. 

The comparison with the observations was performed 
by calculating the x^-statistics for each set of parameters. 
The uncertainty in the derived line intensities due to the 
calibration and the disentanglement of the core and wing 
components was assumed to be 30%. The best fit models 
for the different species agree quite well in the (Tkin, "-H2) 
plane, so all lines are combined into a single estimate 
to constrain the parameters. This is illustrated in Fig. ^| 
where the constraints on densities and temperatures have 
been plotted for given (optimal) values of the column den- 
sities for the individual molecules. For the core component 
an H2 density of 10^ cm~'^ and temperature of 20 K is 
found to be consistent with the observations with a re- 
duced Xrcd of 1.8 for 8 fitted lines. Note that the density 
and temperature are closely coupled, making individual 
determinations somewhat ambiguous, as illustrated in the 

plot, where it is seen that a lower temperature and cor- 
respondingly higher density are equally probable. A lower 
density/higher temperature can also not be completely 
ruled out. For the wing component a best fit density of 
2.0 X 10^ cm~'^ and temperature of 70 K is found with x^cd 
of 3.8 for 13 fitted lines. T he temperature is slightly lower 
than the value quoted bv .Bachiller et all ijlOOS) . but still 
within the mutual uncertainties. The derived column den- 
sities assuming these temperatures and densities are given 
in Table El 

The kinetic temperatures and densities may vary be- 
tween regions traced by different molecular species. Both 
the velocity profiles of the various molecules and the struc- 
ture of their emission in the interferometry maps show 
dissimilarities, indicating chemical differentiation possibly 
due to a combination of the shock evolution and varia- 
tions in the physical conditions. On the other hand as 
illustrated in Fig. [7| the conditions for lines of a particu- 
lar molecule are remarkably homogeneous over the entire 
shock velocity range. It is also found that the derived pa- 
rameters do not vary much when specific molecules are 
included or not. The constraints put on our derived tem- 
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Fig. 10. Best fitted densities and kinetic temperatures de- 
rived from statistical equilibrium calculations with the as- 
sumptions described in the text. In the upper panel the 
results for the core part of the lines are shown and in the 
lower panel results for the wing part. The grey-scaled con- 
tour indicate the la confidence region, while the other con- 
tours, indicate 2a, 3a and 4a confidence levels. For each 
plot the column densities given in Table 13 have been used 
for creating the cut in the cube of models with varying 
temperature, density and column density. For each plot 
the black indicate the best-fit values. The reduced 
is 1.8 (8 fitted lines) for the core component and 3.8 (13 
fitted lines) for the wing component. 



perature and density are in rough agreement with mea- 
surements from other shocked regions from rnolecular out- 
flows Ce.g..lBachiller fc Perez GutierreJlQQTtlCarav et all 
Il998|) . As a first order approximation the two component 
structure therefore seems to describe the excitation condi- 
tions well. Moreover, the derived column densities are well 
constrained with the assumed temperature and density in 
this regime, as illustrated in Fig. El 



Fig. 11. Confidence plots for column density vs. kinetic 
temperature (right) and density (left) for the CS obser- 
vations of the core component of the line profile (upper 
panels) and wing part (lower panels) . The contours corre- 
spond to the la (solid grey), 2a and 3a confidence levels. 
The derived column densities are well constrained for the 
set of temperatures and densities determined from Fig. 1101 
(black '+'). 

Table 3. Column densities for the various molecules from 
statistical equilibrium calculations. 



Molecule 


7V(core)'' 


iV(wing)'' 




[cm-2] 


[cm-2] 


CO" 


5.1x10'' 


1.2x10'' 


CH3OH 


3.4x10'^^ 


5.8x10'^ 


CS 


5.4x10" 


3.2x10" 


HCN 


l.lxlO" 


5.7x10'^ 


HCO+ 


6.6x10'^ 


3.5x10'^ 


H2CO 


8.9x10'^ 


1.8x10" 


N2H+ 


9.0x10'^ 




SiO 




1.4x10'* 


so 


6.2x10" 


1.3x10" 



Notes: "Calculated for = 1 x lO'^cm"^ and T = 20 K (core) 
and T = 80 K (wing). ''Based on '^CO measurements for the 
column densities in the wings and C'*0 measurements for the 
core (assuming a '®0/'*0 isotope ratio of 540). 



5. Discussion 

5.1. Comparison to other protostellar outflows and 
envelopes 

As mentioned in Sect. ^ only a few studies have ad- 
dressed the chemistry in outfiow regions in detail, the 
most thorough being that of the L1157 outflo^w (e.g 
iBachiller &: Perez Gutierred Il997t iBachiller et"all l200l|) 
and t he molecular cond ensation ahead of the HH2 out- 
flows (|Girart et alJ2002|) . These regions differ significantly 
in their context: the CO maps L1157 show an outflow pro- 
gressing through the protostellar envelope and extended 
cloud. This causes a number of bow shocks and conden- 
sations of enhanced density and varying chemistry along 
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the outflow axis. The chemistry in HH2 is thought to 
be induced by irradiation of the molecular condensation 
through UV flux from the bright Herbig-Haro object. In 
such regions especially HCO+ and NH3 should be greatly 
enhan ced according to the models of IViti fc William j 
lll999h. 

Tabled lists the abundances found for the two compo- 
nents of the IRAS2A outflow, calculated as simple ratios 
between the column densities - taking the CO abundance 
relative to H2 at constant value of 1 x 10~^. This means 
that the quoted abundances for the "wing component" are 
averaged over material with a large range of velocities. As 
indicated by the maps and the single-dish line profiles the 
emission from some of the molecules may be more con- 
centrated toward material with lower velocities. In this 
sense, the quoted abundances are therefore lower limits to 
the abundances in the regions where these molecules are 
observed. 

Table ^ also compares the derived abundances 
to those found for the LI 157 and BHR71 outflo\y s 
I Bachiller & Perez Giitierre j Il997t lOarav et all Il998l) . 

the molecular con densation ahead of the HH2 object 
l|Girart et al.ll2002|) and other types of protostellar envi- 
ronm ents - in particular th e IRAS2A protostellar enve- 
lope ij.Targensen et alJl2003bl) the "hot compone nt" of the 
IRAS 16293-2422 envelope ((Schoier et all 120021) and the 
the " C" position of the dark cloud L134N fPic kens et all 
|2000() . The abundances in the IRAS 16293-2422 envelope 
were derived through detailed radiative transfer of the 
dust continuum and molecular line data. The CO, CS and 
HCN abundances quoted are averages over the entire en- 
velope. Since the abundances in the outer part of the en- 
velope may be lowe r due to freeze-out (see for example 
discussion for CO in I J0rgensen et all l)2002|) ). the quoted 
numbers are likely to be lower limits to the abundances in 
the inner, warm region of the envelope. 

In Fig. El the abundances in the two components of 
the IRAS2A outflow, in the IRAS2A envelope and LI 157 
outflow are compared. In particular SiO, CH3OH, SO and 
CS are significantly enhanced by factors 10-10* in the out- 
flow regions compared to the envelope and quiescent dark 
cloud. The abundances in the IRAS2A wing component 
and L1157 agree very well for SO and SiO, but the abun- 
dances of HCN, H2CO and CS are lower in the IRAS2A 
outflow by factors of 10-100. The CH3OH abundances are 
very large for both outflows - between 5 and 10% of that 
of CO. This is especially signiflcant in comparison with 
abundances found in the Orion hot core and the low-mass 
protostar hot core in IRAS 16293-2422. There CH3OH 
is thought to be enhanced through thermal evaporation 
off dust grain ice mantles, but its abundance is a factor 
of 10-100 lower than in the outflow regions. Also the SiO 
abundances are different between the hot core sources and 
the outflow regions, whereas the SO abundances are prac- 
tically identical. For CS and HCN, the abundances have to 
be significantly higher in the inner regions of IRAS 16293- 
2422 in order to match the outflow abundances. Compared 
to the dark cloud, L134N, the CH3OH abundances again 
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Fig. 12. Ab undances in the IRAS2A envelope 
ijjargensen et al.. .2003b) compared to the core and wing 
position components of the I RAS2A outflow (this paper) 
and the L1157 outflow ijBachiller fc Perez Gutierrez! 
I1997I) . The abundances have been normalized to the 
abundances of L134N. For SiO, which has not been 
detected in L134N, a "reference" abundance of 10~^^ has 
been assumed. Note that N2H+ is not detected in the 
outflows, whereas SiO is not detected in the quiescent 
components. 

stand out as significantly increased, together with CS and 
SO. HCN and H2CO have similar abundances in the out- 
fiow region and the dark cloud whereas HCO"'' shows 
slightly lower abundances in the outflow regions. 

5.2. Dynamical time scales 

The dynamical time-scale of the east-west outflow can 
be estimated by assuming that: (1) the terminal veloc- 
ity of the outfiow is equal to the maximal radial velocity 
of the SiO emission, i.e., 29 km s^^, (2) the rest veloc- 
ity of the cloud corresponds to the velocity of the narrow 
non-shocked features (Vlsr — 7 km s~^), (3) the distance 
between the tip of the SiO outfiow and the continuum 
position is the full extent of the outfiow and (4) the in- 
clination between the plane of the sky and the outfiow is 
small, as its high degree of coUimation and extent plus the 
relatively low observed velocities seem to indicate. With 
these assumptions one finds a dynamical age: 



^dyn 



4 X 10^ X ' 



V220pc 



years. 



which is in good agreement with the dynamical timescale 
of 



td 



(3 - 7) X 10^ X 



/ d 



years 



V220pc, 

found bv lBachiller et alJ l)l998() from their CH3OH maps. 

The dynamical timescale calculated this way is sub- 
ject to significant systematic errors. The true extent of 
the outfiow will be Ltruo = iobsCOSi, while the terminal 
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Table 4. Abundances for the two components of the IRAS2A outflow compared to other outflows and protostellar 
environments. 



Molecule 


IRAS2A core" 


IRAS2A wing" 


L1157'' 


HH2'= 


BHR71'' 


IRAS2A env.'= 


IR16293-' 


L134N9 


CO 


=1x10"* 


=1x10-* 


=1x10-" 


=1x10-* 


=1x10-" 


2x10-" 


4x10-" 


=1x10-" 


CH3OH 


6.7x10"'^ 


4.8x10-^ 


2x10-^ 


2x10-* 


2x10-^ 


2x10-3 


3xl0-^t 


8x10-3 


cs 


1.1x10-* 


2.7x10-* 


2x10-^ 


7x10-^° 


6x10-3 


3x10-3 


3x10-3 


1x10-3 


HCN 


2.2x10-3 


4.8x10-3 


5x10-^" 


1x10-3 




2x10-3 


1x10-3 


7x10-3 


HCO+ 


1.3x10-3 


2.9x10-3 


3x10-* 


3x10-* 


9x10-" 


3x10-3 


1x10-3 


8x10-3 


H2CO 


1.7x10-* 


1.5x10-* 


3x10-^ 


2x10-* 




8x10-"' 


6xl0-*t 


2x10-* 


N2H+ 


1.8x10-3 










5x10-3 


IxlQ-" 


6x10-" 


SiO 




1.1x10-'^ 


7x10-* 




7x10-" 


<5xl0-" 


5xl0-3t 


<lxlO-"'' 


so 


1.2x10-* 


1.2x10-^^ 


3x10-^ 


8x10-3 




3x10-3 


3xl0-^t 


6x10-3 



Notes: "-" Molecule not observed. ". . ." Molecule observed but not detected. 



"Core and wing part of IRAS2A outflow (this paper). 
''The L1157 outflow ("Bachiller & Perez Gutierrez 1997") 
■^The HH2 molecular c ondensation ( G i rart et al. 2002) 

''The BHR71 outflow JCarav et al.lll998ll 

"^The IRAS2A protostellar envelo pe (iJorgensen etlolbOOl E)03b^ 

•''The IRAS 16293-2422 envelope dSchoier et alJl2002l) . "f" indicate abundances for the warm inner part of the envelope; other 

abundances are averages over the entire enve lope. 

^The "C" position of the L134N dark cloud jPickens et alJl2000l) 

''Derived using upper limits on [SiO]/[HCN] abundance ratio for L134N of 0.0015 from lZiurvs et alj 



velocity of the outflow relates to the observed maximum 
radial velocity as Vtruc — Kibs sini, modifying the dynam- 
ical timescale by a factor tani. Furthermore the dynami- 
cal timescale at best reflects the properties of the outflow 
at the present moment - changes in the flow velocities 
throughout the history will similarly change the outflow 
dynamical timescale. Thus, it is deflnitely not an unbiased 
indicator of the age of the driving protostellar source. Still, 
it agrees well with the timescales derived from compari- 
son between the envelope structure and collapse models 
( J0rgens en_ct a l. 2003a) and therefore does give an indica- 
tion of the order of magnitude of the appropriate timescale 
to be used when discussing the chemical evolution of the 
shock in the following section. 

5.3. Chemical evolution 

The above analysis shows that a number of species are 
signiflcantly enhanced in the outflow region. The various 
spectral signatures and morphologies in the interferome- 
ter maps together with the varying abundances found in 
the L1157 and IRAS2A outflows indicate different mech- 
anisms regulating the abundances. For the sulfur-bearing 
species, i.e., CS and SO, the enhancement is thought to 
occur as a res ult of enhanced H2S formation as seen in 
hot cores (tPineau des Forets et all.1993: .Charnlev ,1997i|) . 
CS and SO show different morphologies in the interfer- 
ometry maps and different spectral signatures. In partic- 
ular, CS is present in a narrow high velocity component, 
while SO only shows up at lower velocities and with more 
compact emission. In this context the differences and sim- 
ilarities between the abundances found in the LI 157 and 
IRAS2A outflows are also interesting (see Fig. I12() : SiO 
and SO are seen to be quite similar in the L1157 out- 
flow and IRAS2A wing component, while CS, II2CO and 



HCN are found to be 10 to 100 times more abundant in 
the L1157 outflow. An explanation could be differences 
in the abundances of atomic carbon: the production of 
CS is increased with a higher abundance of atomic car- 
bon, whereas the SO is more closely related to the oxygen 
abundance (e.g.. Ivan Dishoeck fc Blakjll998i) . In models 
for gas-phase chemistry, CS, HCN and H2CO are all en- 
hanced with higher abundances of atomic carbon, thus 
potentially explaining the difference between the IRAS2A 
and L1157 outflows. The atomic carbon could either be 
produced in the shock itself or by photodissociation of 
CO in the pre-sho cked gas (e.g., as in the case of IC 443, 
iKeene et al1ll996l) . 

The enhancement of SiO is thought to be caused by 
atomic silicon sputtering from the surfaces of dust grains 
and quickly forming Si O in th e gas-phase through re- 
actions with OH (e.g.. ISchilke et al. 199^ ICaselH et alJ 
ll997HPineau des Forets et al.lll997ll . CH-.OH. H2CO and 
HCN, in contrast, are most likely enhanced through di- 
rect evaporation of ice mantles. Alternative explanations, 
e.g., gas-phase reactions between Cit^ and H2O forming 
CH3OH only produce C H^jOH abundances of 1 - 5 x 10"^ 
(e.g., iMillar et alJll991.) . significantly lower than those 
found in the outflow regions of 10~^ — 10~^ (this paper, 
lRa.chil1er fc Perez Cutierre^ IT997t 
The spectral signatures of CH3OH 



Ra.chi1ler et a,1. 


Il99,^ 


Bachiller et al. 


|1998|) 



compared to SiO indicate that it is formed at lower veloc- 
ities. Furthermore, it is seen that CH3OH peaks slightly 
further downstream compared to, e.g., SiO in the interfer- 
ometer maps. This is clearly illustrated in Fig. 1131 where 
the intensities of SiO, CH3OH and HCO"*" are compared 
along the outflow propagation axis. iGarav et alJ ll200(t 
analyzed regions of different shock velocities in the out- 
flow associated with NGC 2071 and found that CH3OH 
was most prominent in regions with low shock veloci- 
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Fig. 13. Spatial differences between SiO, CH3OH and 
HCO+ along the shock propagation from the interferom- 
eter maps. The emission from each species has been in- 
tegrated in strips perpendicular to an axis aligned in the 
propagation direction of the outflow (position angle of 19° 
with the RA axis and with zero-point at (87" ,-23" ) as in 
Fig. |3J). The SiO and CH3OH data cubes have been re- 
duced to the same spatial resolution as that of the HCO+ 
data. 

ties {vs < 10 km s~^). iGarav et al.l suggested that since 
molecules such as CH3OH, H2CO and HCN are more 
volatile, they would only be capable of surviving at much 
lower shock velocities than would, e.g., SiO. Applied to the 
IRAS2A outflow this explains both the differences in ter- 
minal velocities between SiO on the one hand and HCN, 
H2CO and CH3OH on the other, but also the different 
"onsets" along the propagation direction between CH3OH 
and SiO from the interferometer maps as seen in Fig. 1131 
The characteristic core-wing structure of the observed 
lines is similar to that seen in a SiO survey of protostel- 
lar outflows bv ICodella et al.l l|l999tl . who argued that it 
could be an evolutionary effect with the SiO being pro- 
duced at high velocities and subsequently slowed down 
toward lower velocities. They argued that it would take 
~ 10** years to slow down an outflow-induced shock, 
which is similar to the time it would take SiO to be de- 
stroyed either thro ugh direct accretion onto dust grains 
ijBergin et al . '1998') or through reactions with OH, form- 
ing SiOo |Pineau des Forets et al.l|l997^ . If this picture 
applies to the IRAS2A outflow, it is not surprising that 
SiO has low abundances in the "core component" - or am- 
bient cloud. Since the SiO destruction timescale is similar 
to the dissipation timescale for the protostellar shock, SiO 
is almost completely destroyed in the slow-down phase and 
will therefore not be seen in the low-velocity/quiescent 
component. On the other hand, since SiO is created as a 
direct result of the shock impact l)Pineau des Forets et alJ 
Il997; .Schilke et al, .1997.) , it traces the highest velocities in 
the outflow, together with CO. The characteristic molecu- 
lar depletion timescale at densities of 10^ cm^'^ typical of 
the outflow region (Sect. is on the order of 10"^ years. 



comparable to the outflow dynamical timescale, and could 
therefore explain why, e.g., SiO and CH3OH are not ob- 
served over the entire extent of the outflow back to the 
central protostar. 

The differences between the hot core/warm envelope 
and outflow abundances of, e.g., CH3OH and SiO could 
be caused by differing time scales related to the densities 
in the differing regions: the density in the hot inner part 
of protostellar envelopes is higher by 2-3 orders of magni- 
tude than what is found in the outflow regions. This will 
lead to more rapid destruction of molecules with "anoma- 
lous" abundances, e.g., SiO, either through accretion or 
reactions with other species and therefore also lower abun- 
dances in the envelope regions. Of course the mechanisms 
for producing the given molecules in the first place are 
also likely to be dependent on the environment, further 
complicating the picture. 

The depletion timescale for CH3OH may also be taken 
as an important clock related to the HCO"*" abundance. 
As noted previously, HCO+ stands out compared to the 
other molecules tracing material only in the aftermath of 
the shock. In the L1157 outflow, HCO+ was only found to 
be prominent in the part of the outfl ow close to the driv- 
ing source. Through chemical models. ISergin et"al] l(l998f) 
found that HCO"*" should be destroyed after the passage of 
the shock through reactions with H2O (H2O -I- HCO^ 
H3O"'" -I- CO), but increases later as the water abundance 
reaches lower levels due to freeze-out. This is in fact seen 
in interferometer data as illustrated in Fig. 1131 the emis- 
sion of HCO+ and CH3OH is almost anticorrelated, with 
CH3OH being located closer to the "head" of the outflow 
and HCO+ showing up in the aftermath of the shock. 
As higher abundances of both CH3OH and H2O are ex- 
pected to be results of grain mantle evaporation and the 
timescales for their freeze-out are similar, the HCO"*" and 
CH3OH enhancements should indeed be anticorrelated as 
seen in Fig. E| 

N2H+, like HCO+, is expected to be destroyed by 
reactions with H2O. In contrast to HCO"*", however, 
CO may also be important in destroying N2H+ (e.g., 
iBergin & Langer 1997; Charnley 1991,). Observational 
studies of pre- and proto s tellar objects fe.g..lBergin etall 
l200lt iTafalla et alJliool Ijc^rgensen et al.ll2003a|) suggest 
that N2H+ is enhanced where CO is depleted. The narrow- 
ness of the N2H+ lines and the morphology of the emis- 
sion in the IRAS2A region indicate that this molecule is 
indeed only tracing the ambient cloud material where CO 
may be depleted and not the outflowing material where 
CO is returned to the gas-phase. 

6. Conclusion 

A (sub)millimcter study of the shock associated with 
the NGC1333-IRAS2A outflow has been presented. Both 
single-dish and interferometry line observations are pre- 
sented, which allows for a detailed discussion of both the 
physical and chemical properties in the shocked region and 
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the spatial distribution of emitting species. The main find- 
ings are as follows: 

1. Interferometer observations of the outflow region re- 
veal a distinct morphology with a narrow high veloc- 
ity feature in CS, SiO and HCN, while the low ve- 
locity part traces the more spatially extended mate- 
rial. HCO"'" is seen only in the aftermath of the shock, 
whereas N2H+ shows very narrow lines and does not 
seem to be present in the outflowing gas. 

2. Statistical equilibrium calculations show that the re- 
gion can be divided into two distinct components: low 
velocity material in the ambient, quiescent cloud and 
high velocity material associated with the outflow. It is 
found that the conditions in the quiescent material are 
consistent with a temperature of ~20 K and density 
of 1 X 10® cm^^, while the shocked gas is warmer with 
a temperature of ~70 K and density of 2 x lO^cm"^. 
Within these components, however, the physical condi- 
tions are remarkably homogeneous as indicated by the 
similarities between the lineshapes for different transi- 
tions of various molecules and between the wing emis- 
sion in spectra from interferometry and single-dish ob- 
servations. 

3. The chemistry in the outflow and quiescent regions 
are significantly different. CH3OII, SiO and the sulfur- 
bearing species are significantly enhanced in the out- 
flowing gas, whereas HCN, H2CO, HCO+ and N2H+ 
show abundances more similar to those found in molec- 
ular clouds and protostellar envelopes. Compared to 
the well-studied L1157 outflow the CS, HCN and 
H2CO abundances are markedly lower. This could be 
due to differences in the amount of atomic carbon in 
the shocked or pre-shocked gas. Higher abundances 
of the shock-tracing molecules (in particular SiO and 
CH3OH) in outflows compared to the inner warm en- 
velopes may be related to more rapid destruction of 
these molecules in the envelopes where the densities 
are higher. 

4. A scenario is suggested where the highly collimated 
protostellar outflow is progressing into a region with a 
steep density gradient, focusing the shock and giving 
rise to the narrow morphology observed. This leads to 
a shock-induced chemistry, which can explain both the 
morphologies and qualitatively the abundances of the 
different molecules. In particular, CH3OH is seen to 
be greatly enhanced reaching an abundance of about 
5% of the observed CO abundance in the shocked gas. 
CH3OH is, however, not observed to as high velocities 
as seen in SiO, possibly as a consequence of CH3OH 
being more volatile. Thus together with, e.g., HCN and 
H2CO, CH3OH is not able to survive at higher veloc- 
ities. 

This work illustrates the large impact of protostel- 
lar outflows in shaping the physical and chemical prop- 
erties of their parental environment. The combination 
of high-resolution interferometer observations and single- 
dish spectra makes it possible to address the physical and 



chemical conditions in the shocked and ambient gas and 
to investigate the spatial variation and time-scales charac- 
teristic for the shock induced chemistry. So-far only a few 
shocks have been studied in great chemical detail. Similar 
systematic studies of a large number of different outflows 
will allow for a more detailed comparison between outflows 
and shocks of different velocities and energetics and in 
different environments. Future observations with facilities 
such as the SMA, CARMA, and ALMA will allow further 
studies of the variation of physical and chemical conditions 
in shocks through high resolution, high sensitivity multi- 
transition molecular line observations. Also high spatial 
resolution observations of H2O lines with Herschel-HIFI 
can conflrm the anticorrelation between HCO+ and H2O. 
All such more detailed observational studies will serve as 
important starting points for more detailed physical and 
chemical models for shocks in protostellar environments. 
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